Coherence and Raman sideband cooling of a single atom in an optical tweezer 
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We investigate quantum control of a single atom in an optical tweezer trap created by a tightly 
focused optical beam. We show that longitudinal polarization components in the dipole trap arising 
from the breakdown of the paraxial approximation give rise to significant internal-state decoherence. 
We show that this effect can be mitigated by appropriate choice of magnetic bias field, enabling 
Raman sideband cooling of a single atom close to its three-dimensional ground state in an optical 
trap with a beam waist as small as w = 900 nm. We achieve vibrational occupation numbers of 
n r = 0.01 and n a = 8 in the radial and axial directions of the trap, corresponding to an rms size 
of the atomic wavepacket of 24 nm and 270 nm, respectively. This represents a promising starting 
point for future hybrid quantum systems where atoms are placed in close proximity to surfaces. 



Single atoms in "optical tweezer" traps [1 are a 
promising resource for various applications in quantum 
science and engineering. They can be individually moved 
[2] , manipulated [3j [4] , and read-out [5j in a manner sim- 
ilar to trapped ions. At the same time, they may be 
strongly coupled to photonic [SI [7], plasmonic [8j, or other 
solid-state systems [9HTT]. opening a new frontier for the 
realization of quantum networks and hybrid quantum 
systems. These intriguing applications require trapping 
single ultra-cold atoms near surfaces at distances well be- 
low an optical wavelength. While this is challenging for 
ions [12 , and magnetically trapped atoms [H[T3], it is 
readily achievable with neutral atoms in optical dipole 
traps. 

The collisional blockade regime pQ of optical dipole 
traps is an attractive starting point for such experiments 
because it provides a simple way to load and tightly con- 
fine single atoms starting with only a conventional optical 
molasses. In several experiments, it has been found that 
the temperature of an atom loaded from a molasses into 
an optical dipole trap in the collisional blockade regime 
is in the range 30 fiK to 180 /iK [3 H CD UMI]. This 
elevated temperature compared to free-space cooling has 
been identified as a limiting factor in many recent ex- 
periments, as the thermal motion reduces the coherence 
time [3j [4j [16] and impedes full quantum control [T7J [18] . 
Moreover, lower temperatures and a reduction in the spa- 
tial extent of the atomic wavepacket are necessary to con- 
trol the atom at sub- wavelength distances from a surface, 
or to implement proposed quantum gates using collisional 
interactions between two ground state atoms [T9] . 

One major challenge to laser cooling and quantum con- 
trol are polarization effects associated with the break- 
down of the paraxial approximation in very tightly fo- 
cused optical dipole traps. Vector light shifts arising 
from elliptical light polarization [20] are known to be a 
major obstacle to cooling and manipulating atoms in op- 
tical dipole traps. In the paraxial limit, the vector light 
shift can be eliminated by using a linearly polarized trap- 



ping beam. However, in the tightly focused regime, non- 
paraxial effects produce a longitudinal and spatially in- 
homogeneous polarization component near the focus that 
cannot be eliminated. The resulting state-dependent 
trapping potential leads to dephasing of internal state 
superpositions [21] and fluctuating dipole force heating, 
which impairs internal state manipulation as well as laser 
cooling. 

In this Letter, we present a detailed study of the lon- 
gitudinal polarization component of dipole trap formed 
by a high numerical aperture lens, demonstrate how the 
undesirable effects arising from the longitudinal polar- 
ization can be partially compensated using a properly 
oriented magnetic bias field, and apply these results to 
perform Raman sideband cooling of a single atom. After 
cooling, the atom is in the ground state along the two ra- 
dial directions (n r = O.OIIq oi)' anc ^ occupies just a few 
quantum states (n a = 8.1(8)) in the axial trap direction. 
The corresponding rms size of the atomic wavepacket is 
given by the ground state length of 24 nm in the radial 
directions, and a thermal extent of 270 nm in the axial 
direction, which represents a hundred-fold reduction in 
spatial volume, and a reduction by 10 4 in phase-space 
volume, over the starting conditions. 

The origin and effects of the longitudinal polarization 
component can be understood in the framework of ray 
optics (see Figure []Ji). The light entering a lens con- 
sists of parallel rays with linear polarization transverse 
to their propagation direction. Upon passing through the 
lens, rays are refracted according to their distance from 
the optical axis, and their polarization directions must 
also deflect to remain transverse to the ray [22]. In the 
diffraction-limited volume around the focus, all of these 
rays interfere and the resulting field is elliptically polar- 
ized. Following Fig. [T^i, two features of the polarization 
near the focus emerge: (1) the polarization vector is ro- 
tating in the plane set by the incident polarization vector 
and the optical axis, and (2) the sense of this rotation is 
opposite above and below the optical axis. 
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FIG. 1. (a) Diagram showing the origin of elliptical polarization near the focus (see text), (b) Cut through the focal plane 
for a — 0.43. Contour lines show C x , which is C x scaled to the local intensity \E{f)\ 2 /\E{fmax)\ 2 • Shading shows Gaussian 
intensity profile, (c) Dephasing rate between the states |1) and |2) as a function of bias field along the Z axis, with a trap 
wavelength A = 815 nm. The improvement at large bias fields is due to suppression of the polarization gradient. Fit is to 
model described in text: 770 + r\ c %rc are background dephasing rates from the finite detuning and slight elliptical polarization of 
the dipole trap; rj pg arises from the longitudinal polarization. Inset: Ramsey measurement of dephasing rate between |1) and 
|2) at B 10.5 G. 



For light that is far-detuned compared to the excited- 
state hyperfine structure, the vector light shift for alkali 
atoms in the ground state is [20j [23] : 

U(T) = -U (r)p^-C(r)-g F F (1) 

where Uq(t) is the scalar dipole trap potential, Si and 
#2 are the detunings from the Dl and D2 lines, respec- 
tively, e(r) is the local (unit norm) polarization vec- 
tor, F is the total angular momentum operator and 
Qf = [F(F + 1) - 1(1 + 1) + J (J + 1)] /F(F + 1). The 
term C = Im [e(r) x e*(r)] is a basis-independent way of 
expressing the normal vector of the polarization ellipse 
and the degree of ellipticity (with magnitude ±1 for cir- 
cularly polarized light; for linear polarization). Using 
the vector Debye integral [22] , we have numerically com- 
puted the polarization near the dipole trap focus (Fig. 
[Tj}). The most important term is the polarization gradi- 
ent dC x /dy. For a lens with numerical aperture a, the 
maximum gradient, occurring at the beam focus, is well 
approximated by 3.1a sin a/A if the lens aperture is uni- 
formly illuminated, and 2.6a sin a/ A if the illumination 
is a Gaussian beam with a 1/e 2 diameter equal to the 
lens aperture diameter. 

In the experiments presented here, we use an optical 
dipole trap for 87 Rb atoms operating at At = 815 nm 
with a depth U$ = 0.82 mK. Independent measurements 
of the depth and radial trap frequency (uj r = 2tt x 100 
kHz) allow us to extract a 1/e 2 beam radius w = 900 nm 
and an effective lens aperture a = 0.43. For these pa- 
rameters, we expect dC x /dy = 0.57/ am. Since the state- 
dependent potential in Equation (|l| is linear in F, it 
produces the same energy shifts as a magnetic field, and 
dC x jdy can also be expressed as an effective magnetic 
field gradient with magnitude B' x = 1.4 G/jiva at the 



trap center. This gradient gives the trapping potential a 
significant state-dependent component. 

Specifically, in the absence of an externally applied 
magnetic bias field, atoms in different magnetic sub- 
levels experience trapping potentials that are displaced 
by Ax = iiB^(gF^F)B' x / '(racj 2 ), where /ieA^^rnp) is 
the difference in the magnetic moment between the two 
sublevels. For A(^t^f) = 1/2, Ax = 11 nm, which 
is not negligible compared to the ground state length 
^h/2muj = 24 nm. While this state-dependent displace- 
ment could be useful for Raman cooling or other motional 
state manipulations along this axis [24j [25] , it also leads 
to rapid internal-state decoherence on the timescale of 
the radial trap oscillation period. 

To circumvent this problem, we can apply a bias mag- 
netic field in a direction orthogonal to x. In this case, 
the effective field gradient is suppressed since B tot = 
yjBl + (B f x y) 2 ~ B + (B /2 /2B )y 2 , and the gradient 
results in only a state-dependent change in the strength 
of the harmonic trap potential. Quantitatively, super- 
positions of magnetic sublevels that experience different 
trapping potentials of the form Ui(r) = (1 + rf)U2 (r) are 
dephased with a coherence time XJ = 0.97 x 2h/(kBTr]) 
[21], where T is the temperature of the atom and fc# is 
the Boltzmann constant. In the presence of a large or- 
thogonal bias field, the polarization gradient contributes 
to r] as rj P g = fiBA(gFmF)B x /(3muj 2 Bo) (the factor of 
1/3 results from averaging over the three trap axes). We 
can use the dependence on Bq to accurately measure B\ 
and also to suppress the dephasing with large Bo. 

We measure the decoherence between the states |1) = 
\F = l,ra F = -1) and |2) \F = 2,m F = -2) using a 
sequence that starts with loading a single atom into a 
tweezer trap with a depth of 1.6 mK at zero bias field, 
then ramping down the trap depth to 0.82 mK as we 
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ramp up the bias field Bq to the desired value. The atom 
is optically pumped into |2), the hyperfme transition is 
driven by a two-photon Raman process in a Doppler-free 
configuration, and the state detection is accomplished 
using a push-out beam; these details are described in 
more detail below. T 2 * is extracted from a Ramsey-type 
measurement, using a fit to the function introduced in 

ED- 

At each trap wavelength, we fit 1/T 2 * = 
1.03 (770 + Varc + Vpg) (kBT/2H). The only free parame- 
ters in this fit are the degree of circular polarization in 
the incident dipole trap beam due to birefringence (r] C i rc ) 
and the strength of the effective field gradient B' x . The 
temperature is determined independently (T = 40/iK 
for this measurement, see below for technique). 770 
reflects the different trapping potentials for F = 1 and 
F = 2 atoms due to the finite trap detuning, and is 
approximately given by the ratio of ground state hyper- 
fine splitting to trap detuning. At a trap wavelength 
of (802,815) nm, we find B' x = (2.4,1.4) G//im, corre- 
sponding to gradients dC x /dy = (0.46(6), 0.54(3) )//im, 
in reasonable agreement with our estimate of 0.57 / jira. 

Having developed a detailed understanding of trap- 
induced decoherence in this system, we now turn to Ra- 
man sideband cooling. We use three orthogonal running- 
wave fields to drive Raman transitions (R1-R3, see Fig. 
^l), as in trapped ion experiments [26]; the different fre- 
quencies are generated by means of an electro-optic mod- 
ulator [27]. Optical pumping to the |2) state is provided 
by circularly polarized beams propagating along the mag- 
netic field axis, addressing the F = 1 — > F' = 2 and 
F = 2 F' = 2 transitions on the D2 line. The fre- 
quencies of the lasers are set to the measured resonances 
in the dipole trap, which are shifted by ~ 30 MHz from 
the resonances in free space. The intensities of the two 
beams about 100 times less than saturation. This ensures 
that the atom only scatters photons elastically, so heat- 




FIG. 2. (a) Schematic of optical paths. The dipole trap is 
formed at the focus of an aspheric lens (AL), with a 1/e 2 
radius of 900 nm. Beams Rl, R2 (both circularly polarized) 
and R3 (linearly polarized orthogonally to Bo) drive Raman 
transitions. The detuning of the Raman beams A = 13 GHz. 
Optical pumping is along the path marked OP (circularly 
polarized), (b) States in 87 Rb used for Raman transitions 
and optical pumping. 
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FIG. 3. (a) Release and recapture temperature measure- 
ment. (Closed, open) circles show measurements (before, af- 
ter) radial cooling. A Monte Carlo model yields kinetic en- 
ergies K such that 2K/kB — 52(4) /iK before cooling, and 
(2K r /k B ,2K a /k B ) = (2.4(1), 158(14))/iK after cooling. (b,c) 
Measurement of the axial kinetic energy before and after cool- 
ing the axial mode, by Doppler spectroscopy, (b) After radial 
cooling only, 2K a /kB = 129(19)/xK. (c) After radial and axial 
cooling, 2K a /k B = 8.1(l)/xK. 



ing due to fluctuating dipole forces associated with the 
(anti-) trapping potential for the excited state is avoided 
[28] . For diagnostic purposes, we measure the F = 1 pop- 
ulation by pushing out the atoms in F = 2 using a beam 
resonant with the F = 2 — >> F r = 3 transition on the 
D2 line, then measuring whether the atom has remained 
trapped by turning the molasses back on. This beam is 
circularly polarized and propagates along the same path 
as the optical pumping beams. 

In a typical experiment, we load an atom from the 
MOT into the optical dipole trap with a depth of 1.6 
mK at zero bias field, then decrease the trap depth to 
0.82 mK while ramping the bias field Bq up to 7.5 G. 
Lowering the trap depth serves to increase the coherence 
time while leaving the trap frequencies high enough that 
sideband cooling is still achievable, with (uj ri uo a ) = 2tt x 
(100, 15.6) kHz. All temperatures reported in this paper 
are measured in the 0.82 mK deep trap. We cool the 
atoms in the following sequence: we first apply the R2 
and R3 beams (Fig. [2| and the optical pumping beams 
together for 10 ms to continuously cool the radial modes; 
then, we perform ten cycles consisting of 2 ms of axial 
cooling using the Rl and R2 beams, followed by 4 ms of 
radial cooling using the R2 and R3 beams again. This 
sequence prevents the radial modes from heating while 
the axial cooling proceeds. 
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FIG. 4. Sidebands showing final temperature in the (a) radial 
and (b) axial directions. In (a), the red and blue sideband am- 
plitudes are fit to independent lorentzians; their ratio yields a 
radial temperature n = O.Ol^o'.oi- I nse ^: same measurement 
with shorter pulse length so the carrier is also resolved. In 
(b) 9 peaks are fit with independent heights, but equal spac- 
ings and widths. The heights are well-described by a thermal 
distribution with n a — 8.1(1). 



The parameters for the first radial cooling phase are 
optimized by measuring the temperature using a release 
and recapture technique [29]. This data, shown in Fig- 
ure [3^i, is fit using a Monte-Carlo simulation [14J. The 
initial kinetic energy K is such that 2K/kB = 52/iK; 
the measurement after cooling yields anisotropic kinetic 
energies of 2K r /kB = 2.4(1 )/iK in the radial direction 
and 2K a /kB = 158(14)/iK in the axial direction (the 
release and recapture technique is only weakly sensitive 
to the axial mode). The fitted kinetic energies repre- 
sent the global minimum in y 2 over the entire space of 
three independent energies for each axis, including un- 
physical temperatures less than the ground state energy 
ftjjj/2kB = 2A\±K for the radial modes. The agreement 
of the measured kinetic energy with that of the zero- 
point motion suggests that we have reached the radial 
ground state after this cooling phase alone. The radial 
cooling works best with a two-photon Rabi frequency 
&R2,R3 = 2tt x 17 kHz and a detuning of —2tt x 100 kHz 
from the two-photon resonance, corresponding to — uo ri 
as expected. 

To characterize the axial temperature independently 
after the radial cooling, we measure the Doppler width 



of the |2) to |1) transition when driven with the Rl and 
R2 beams. The wavevector Aki 2 = k^i — k# 2 has a 
projection onto the axial and radial directions, but the 
Doppler profile should mostly be sensitive to the axial 
mode here since the radial degrees of freedom are already 
cold. After the first stage of radial cooling, we measure a 
kinetic energy of 2K a /k B = 129(19)/xK (Fig. |3]d). After 
optimization, we obtain a feature with a width corre- 
sponding to 2K a /k B = 8.1(l)/iK (Fig. This data 
is fitted to a Gaussian, which conservatively assumes no 
power broadening. The optimum cooling parameters are 
a two-photon Rabi frequency of ^#1,^2 ~ 27r x 5 kHz 
and a detuning of — 2tt x 60 kHz. The parameters used 
for the interleaved radial cooling phases are the same as 
above. 

To obtain more precise measurements of the final tem- 
perature of the atom, we resolve the asymmetric motional 
sidebands along two axes. The ratio of the sideband am- 
plitude gives information about the vibrational state oc- 
cupation of the atom [26] . 

Figure [4^i shows the sidebands measured in the ra- 
dial direction with small ^#2,^3 • The blue sideband 
is essentially absent, with a fitted amplitude 100 times 
smaller than the red sideband. From this, we extract 
a final temperature for the radial degrees of freedom of 
n = 0. Ol^o'oi- We do not know to what extent the two 
radial modes are non-degenerate or what the preferred 
axes are, but from the release- and-recapture data show- 
ing that both modes must be very cold, and the fact that 
the spectrum shown here does not change if we measure it 
at a different time after the cooling (up to 100 ms later), 
we infer that the two modes are not perfectly degenerate 
and the R2+R3 beams address both modes. Therefore, 
we conclude that this spectrum reflects the temperature 
of both radial modes. 

We also resolve the axial motional sidebands using 
the Rl and R2 beams at very low power, and observe 
a spectrum with nine peaks that is slightly asymmet- 
ric (Fig. [4)3). We find that the ratios of the measured 
peak heights correspond very well to a thermal distribu- 
tion p nn oc exp(— n/n a ) with a mean vibrational number 
n a = 8.1(1). The corresponding energy (n a + l/2)hu a = 
6.5 jaK xks is similar to the result of the Doppler mea- 
surement above. 

Several properties of the cooled atom are worth noting. 
The heating rate for the radial degrees of freedom is very 
low: we observe heating An < 0.3 over 200 ms, which is 
consistent with what is expected from photon scattering 
only. We also do not observe radial heating if we translate 
the atom over distances ~ 20/im in ~ 10 ms using a scan- 
ning galvanometer mirror. Lastly, the fraction of atoms 
lost from our dipole trap (~ 30 %) during the Raman 
cooling time (150 ms) is due entirely to background gas 
collisions (P ~ 10 -8 torr), and not to the cooling process 
itself. We find that decreasing the Rabi frequency Qri^ 
and detuning during the last cooling phase does not de- 
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crease the final axial temperature. This is possibly due 
to one of the following considerations: (1) in our appara- 
tus, the R1+R2 beam pair couples motion along different 
axes, making it difficult to isolate the axial direction; (2) 
the ground state Lamb-Dicke factor is close to one on this 
axis (kAz = 0.72); and (3) the optical pumping beams 
are oriented along the axial direction, potentially causing 
extra heating which we have not characterized. At the 
same time, we are not aware of any fundamental effects 
that would prevent cooling to the ground state in this 
system. 

We have quantified the coherence properties of single 
atoms subject to the polarization distortions present in 
tightly-focused and near-field optical dipole traps, and 
shown that the unwanted effects may be mitigated by ap- 
plying a magnetic bias field in the appropriate direction. 
Using this result, we performed Raman sideband cooling 
of a single atom in an optical tweezer trap close to its 
ground state, achieving final mean occupation numbers 
(n r ,n r ,n a ) = (0.01,0.01,8.1). This technique should 
find immediate application to a variety of experiments 
using optical dipole traps. 

We acknowledge funding from the NSF, CUA, 
DARPA, AFOSR, MURI and the Packard Foundation. 
JDT acknowledges support from the Fannie and John 
Hertz Foundation and the NSF GRFP. 

Note After completion of this work, we have become 
aware of a related demonstration of Raman sideband 
cooling in an optical tweezer [30] . 
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